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This  paper  discusses  the  testing  of  a  single-phase  spray  cooling  system  that  was  flown  on 
the  NASA  KC-135  Reduced-Gravity  Research  Aircraft.  An  experimental  package, 
consisting  of  a  spray  chamber  coupled  to  a  fluid  delivery  loop  system,  was  fabricated  for 
variable  gravity  flight  tests.  The  spray  chamber  contains  two  opposing  nozzles  spraying  on 
Indium  Tin  Oxide  (ITO)  heaters.  These  heaters  are  mounted  on  glass  posts,  which  are  part 
of  a  sump  system  to  remove  unconstrained  liquid  from  the  test  chamber.  Thermocouples 
mounted  in  and  around  the  posts  were  used  to  determine  both  the  heat  loss  through  the 
underside  of  the  ITO  heater  and  the  heat  extracted  by  the  spray.  During  flight  tests,  for 
Weber  numbers  of  We  =  771  ±  19  and  757  ±  15,  the  non-dimensional  heat  input  was  varied 
from  GA  =  30  to  110  for  the  non-dimensional  grouping  (Fr1/2Ga)12  =  20  to  66.  Flight  test 
data  and  terrestrial  data  were  compared  to  analytical  and  numerical  solutions  in  order  to 
evaluate  the  heat  transfer  in  the  heater  and  support  structure.  In  general,  the  Nusselt 
number  at  the  heater  surface  was  found  to  decrease  with  increasing  (FrI/2Ga)12. 


a 

b 

Bi 

A, 

/ 

Fr 

g 

G 

GA 

Ga 


Nomenclature 

acceleration  level,  g 

radius  of  the  glass  heater  post  assembly,  m 

Biot  number,  hmn  b/kpost 

droplet  diameter,  m 

heater  conduction  loss  fraction 

Froude  number,  v2/aDp 

volumetric  heat  generation,  W/m3,  Q/n/rS 

gb2/(Too  ,wall  ^liti) 

non-dimensional  heat  input,  Q/( n  b  £htr  TKjWan) 

Galileo  number,  aZ)p3p2/p2 

spray  heat  transfer  coefficient  at  the  top  of  the  heater,  W/(m2-K) 
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A  wail  =  heat  transfer  coefficient  along  the  side  of  the  post,  W/(m2-K) 

H  =  glass  layer  height,  m 

Hp  =  height  of  the  glass  post  assembly  and  heater,  m 

Jo,  J\  =  Bessel  functions 

k  =  thermal  conductivity,  W/(m-K) 

K  ^fluic/^htr 

L  =  height  of  the  glass  post  assembly,  m 

Nu  =  Nusselt  number,  htop  b/knuid 

q  =  heat  flux,  W/m2 

q  i  =  heat  flux  lost  due  to  conduction  into  the  glass  post  assembly,  W/m 2,fQ/nb2 , 

Q  =  heater  input  power,  W 

r  =  radial  coordinate,  m 

R  =  non-dimensional  radial  coordinate,  r/b 

t  =  time,  s 

T  =  temperature,  °C 

Tint  =  interface  temperature,  TCI,  °C  (located  at  z  ^  / A  +  //4  +  II J 

T,  =  heater  surface  temperature,  °C 

^sat  =  saturation  temperature,  °C 

T^op  ~  local  freestream  temperature  of  liquid  flowing  over  heater  surface,  °C 

T’cB.waii  =  local  freestream  temperature  of  liquid  flowing  along  post  wall,  °C 

v  =  droplet  velocity,  m/s 

V  =  volumetric  flow  rate,  m3/s 


n  =  index  integers 

post  =  region  2,  glass  post  assembly 

uf  =  upper  film 

2, 3, 4, 5  =  layers  2-5  of  the  glass  post  assembly 
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I.  Introduction 


THE  dissipation  of  thermal  energy  from  current  and  future  electronics  for  proposed  high  power  space  and  air 
vehicle  applications  has  resulted  in  the  development  of  thermal  management  approaches  tailored  for  high  heat 
flux  acquisition  and  high  thermal  energy  transport.  Compounding  the  already  difficult  thermal  problem  is  the  fact 
that  these  heat  generating  sources  generally  have  low  thermal  masses.  This  requires  special  attention  be  given  to  the 
development  of  robust  thermal  management  approaches  that  must  operate  reliably  in  space  and  aircraft  applications 
while  being  subjected  to  variable  body  forces  resulting  from  a  variable-gravity  environment.  Recently,  the  emphasis 
has  been  on  the  development  of  two-phase  cooling  concepts  to  minimize  system  penalties  as  compared  to  the 
simpler  single-phase  cooling  approaches.  Several  two-phase  thermal  management  approaches  have  been  or  are 
being  considered  for  space  and  air  vehicle  environments.  These  include  two-phase  flow  and  heat  transfer1'2,  pool 
boiling3'4'5,  and  spray  cooling6'7. 

Many  variable  gravity  studies  have  focused  on  the  physics  involved  in  two-phase  flow  with  heat  transfer,  and 
experiments  have  been  performed  in  drop  towers,  reduced-gravity  aircraft  flights,  and  on-orbit.  Two-phase,  one- 
component  flow  with  heat  transfer  in  microgravity  is  seen  in  many  thermal  management  systems  such  as  flow 
boiling  systems,  heat  pipes,  loop  heat  pipes,  and  capillary  pumped  loops.  Work  has  been  done  to  try  to  use 
similarity  considerations  to  scale  terrestrial  models  to  predict  the  behavior  of,  and  to  design,  prototypes  for 
microgravity1'2. 

Variable  gravity  research  on  pool  boiling  with  and  without  subcooling,  concentrating  on  microgravity 
conditions,  has  been  performed  by  many  researchers3'4'5.  It  has  generally  been  found  that  microgravity  conditions  do 
not  significantly  change  the  heat  transfer  coefficient  for  nucleate  boiling.  However,  the  Critical  Heat  Flux  (CHF)  is 
generally  reduced  in  microgravity.  This  reduction  in  the  CHF  is,  in  part,  a  result  of  the  variation  in  bubble  dynamics 
in  microgravity  due  to  the  lack  of  buoyancy  forces.5  Observations  by  Kim  et  al.4  showed  that  small  bubbles  will 
coalesce  into  a  large  bubble  on  the  surface  of  the  heater  in  reduced  gravity.  During  subcooled  boiling, 
thermocapillary  flows  can  develop  around  bubbles  forcing  the  bubble  to  the  heater  surface  while  drawing  warm 
liquid  away  from  the  heater  and  supplying  cool  liquid  to  the  heater/ 

Many  researchers  have  investigated  two-phase  sprays  for  the  thermal  management  of  devices  generating  high 
heat  fluxes;  however,  there  has  been  little  research  addressing  the  physics  and  ultimate  performance  of  two-phase 
spray  cooling  in  variable  gravity  environments.  The  development  of  concepts  such  as  spray  cooling  in  these 
environments  must  be  supported  with  a  sound  understanding  of  the  micro-  and  macro-scale  two-phase  spray 
thermophysics. 

Several  variable-gravity  spray  cooling  experiments  were  performed  by  Yoshida  et  al.6  and  Kato  et  al.7  onboard  a 
reduced-gravity  aircraft.  The  experiments  principally  looked  at  transient  and  steady-state  spray  cooling  in  the 
microgravity  portion  of  the  reduced  gravity  aircraft  flights.  A  copper  block  heater  was  used  for  transient  heat 
transfer  experiments,  and  a  clear  thin  fdm  ITO  heater  on  a  clear  base  was  used  for  steady-state  experiments  with 
flow  visualization.  The  volumetric  spray  fluxes  used  were  between  F"  =  1.36  x  10~4  to  3.70  xlCT4  m3/m2-s.  These 
fluxes  were  low  enough  that  a  continuous  film  of  coolant  was  not  allowed  to  develop  on  the  heater  surface.  The 
Weber  numbers  for  their  experiments  ranged  from  We  =  28  to  622,  with  all  but  two  of  the  experimental  cases  having 
Weber  numbers  below  We  =  60. 6  Fluid  management  issues  that  would  be  encountered  for  a  continuous-flow, 
closed-loop  microgravity  spray  cooling  system  were  not  addressed.  No  attempt  was  made  during  the  experiments  to 
retrieve  and  reuse  the  working  fluid  in  a  closed-loop  system. 

Baysinger  et  al.8  presented  the  design  and  preliminary  testing  of  a  variable-gravity  spray  cooling  experiment. 
To  build  on  the  work  of  the  previously  mentioned  researchers,  the  primary  purpose  of  the  study  by  Baysinger  et  al. 
was  to  investigate  heat  transfer  and  fluid  management  issues  for  a  continuous-flow,  closed-loop  spray  cooling 
system  subjected  to  a  variable-gravity  environment.  Tests  were  conducted  onboard  the  NASA  KC-135  Reduced- 
Gravity  Research  Aircraft,  which  provided  the  variable-gravity  environment  by  following  a  parabolic  flight 
trajectory.  Emphasis  was  placed  on  developing  an  initial  transient  analytical  model  to  predict  temperatures  and  the 
spray  heat  transfer  coefficient  within  the  heated  region  of  the  experiment  for  low  heat  loads.  The  flow  rate  was 
sufficient  to  keep  the  heater  surface  completely  immersed  in  the  working  fluid. 

The  objective  of  the  present  investigation  is  to  continue  the  work  presented  by  Baysinger  et  al.8,  and  to 
investigate  the  effects  of  variable-gravity  on  single-phase  spray  cooling  by  varying  the  non-dimensional  heat  input 
(GA  =  30  to  110)  for  the  non-dimensional  grouping  (Frl  2Ga)12  =  20  to  66.  Fluid  management  and  heat  transfer 
thermophysics  were  investigated  for  Weber  numbers  We  =  771  and  We  =  757. 
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II.  Experimental  Design 


The  experimental  apparatus  consists  of  two  primary  components  similar  to  that  described  by  Baysinger  et  al.s 
The  first  component  was  a  spray  test  chamber  containing  two  glass  heater  posts,  two  spray  nozzles,  a  primary 
condenser  surface,  and  two  sump  configurations  to  collect  the  liquid  and  condensate.  The  second  component 
consisted  of  a  flow  loop  to  manage  the  working  fluid  (FC-72)  and  a  water  loop  used  to  control  the  temperature  of  the 
spray  test  chamber.  A  photograph  of  the  overall  experimental  apparatus  is  shown  in  Fig.  1(a).  Figure  1  (b)  shows  the 
spray  test  chamber,  which  was  a  cylinder  having  interior  dimensions  20.2  cm  (7.97  in)  in  diameter  and  15.2  cm  (6 
in)  in  length.  It  was  made  of  standard  vacuum  components  to  provide  a  leak-tight  environment  for  the  spray-cooling 
experiment. 

Each  glass  heater  post  assembly  was  cylindrical  with  a  diameter  of  0.016  m,  and  consisted  of  five  layers,  as 
shown  in  Fig.  2.  The  first  layer  was  a  thin-film,  transparent,  indium-tin  oxide  (ITO)  heater,  which  was  sputter- 
deposited  onto  the  second  layer,  a  0.5-mm-thick  fused  silica  substrate.  This  combination  was  chosen  due  to  its  low 
thermal  mass  and  to  match  the  coefficient  of  thermal  expansion  as  closely  as  possible  between  the  ITO  layer  and  the 
glass  layer  to  prevent  delamination.  The  low  thermal  mass  of  the  ITO  heater  resulted  in  a  rapid  response  time  such 
that  the  system  could  reach  steady-state  conditions  during  the  parabolic  flight  profile,  which  typically  lasted  from  20 
to  25  s.  The  heater/ substrate  was  mounted  onto  a  1 -mm- thick  wafer  of  borosilicate  glass  using  optical  cement.  This 
subassembly  was  then  adhered  to  an  identical  1-mm-thick  wafer  of  borosilicate  glass,  which  was  finally  cemented  to 
the  main  glass  post,  which  was  also  made  of  borosilicate  glass.  The  purpose  of  the  two  intermediate  layers  was  to 
allow  for  channels  to  position  thermocouples,  as  shown  in  Fig.  2(b).  The  heater  post  assembly  also  allowed  for  an 
optical  path  for  visual  observation  of  the  heater  surface.  A  preliminary  analysis  was  used  to  predict  temperatures 
throughout  the  post  in  order  to  determine  the  best  placement  for  the  thermocouples.  The  glass  heater  post  assembly 
was  integrated  into  a  sump  configuration,  located  adjacent  to  each  nozzle,  and  served  to  remove  excess  fluid  from 
the  chamber  and  spray  surface  through  an  annulus  between  the  heater  post  and  sump  wall,  as  shown  in  Fig.  3. 


b) 


Figure  1.  Experimental  flight 
package  showing  the:  (a) 
complete  package  as  flown  and 
( b )  test  chamber. 


b) 


Figure  2.  Magnified  view  of  the 
heater  region  showing  the;  (a) 
side  view  with  the  glass  layers 
and  (b)  top  view  with 
thermocouples. 


b) 


Figure  3.  Assembled  glass  heater 
post  showing;  (a)  thermocouples 
installed  and  ( b )  a  Solidworks 
drawing  of  the  assembled  sump 
with  the  glass  heater  post. 


Table  1.  Nozzle  performance  characteristics. 


Volumetric  Flow  Rate,  V 
(m3/s)/(gph) 

Nozzle  Pressure 
(kPa)/ (psig) 

Droplet  Diameter,  Z)p 
(pm) 

Droplet  Velocity,  v 
(m/s) 

7.36><10'6  /  7.0 

248/36 

48 

9.0 

l.OlxlO'5  /  9.6 

324/47 

48 

12.0 

1.21x10  s  / 11.5 

469/68 

47 

15.0 
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The  spray  nozzle  characteristics  were 
evaluated  using  a  two-axis  phase  Doppler 
anemometer  to  determine  spray  uniformity  and 
droplet  size  distribution  as  described  by 
Baysinger  et  al.s  Typical  droplet  velocities, 
mean  droplet  diameters,  and  nozzle  pressures 
are  shown  in  Table  1.  The  distance  from  the 
nozzle  to  the  heater  surface  was  set  to 
approximately  1  cm.  This  distance  was  based 
on  a  spray  cone  angle  of  60°  such  that  the 
spray  contacted  only  the  heater  surface  with 
minimal  overspray.  The  design  and  placement 
of  the  nozzle  in  relation  to  the  heater  cap, 
heater,  and  annular  sump  entrance,  shown  in 
Fig.  3{b),  allowed  for  vapor  generated  at  the 
heater  surface  to  be  transported  past  the  nozzle 
into  the  test  chamber  interior  or  entrained  with 
the  liquid  through  the  sump  annulus.  In 
addition,  excess  liquid,  either  due  to  droplet 
bounce  or  flooding,  could  be  transported  past 
the  nozzle  into  the  test  chamber  interior. 

Excess  vapor  and  liquid  in  the  chamber 
were  controlled  by  an  annular  wick  structure, 
both  to  aid  in  condensation  and  to  contain  the 
condensate,  as  shown  in  Fig.  1  (b).  The  wick 
structure  was  formed  by  lining  the  internal 
chamber  wall  with  one  layer  each  of  100  and 
150  stainless  steel  mesh,  with  the  150  mesh 
being  nearest  the  wall.  An  additional 
cylindrical  screen  was  placed  in  the  chamber 
interior  to  provide  for  an  annular  space 
between  the  screen  and  the  chamber  wall.  This 
annular  space  between  the  chamber  wall  and 
inner  screen  was  filled  with  stainless  steel  wool  to  help  confine  the  unconstrained  liquid  in  the  microgravity 
environment.  In  addition,  a  pressure  transducer  was  installed  to  measure  the  vapor  pressure  within  the  chamber  but 
external  to  the  annular  wick  structure.  Transient  variation  in  acceleration  throughout  the  flight  profile  was 
monitored  using  a  tri-axis  accelerometer,  with  an  uncertainty  of  ±0.03-g,  mounted  onto  the  test  chamber.  For  this 
investigation,  the  lower  nozzle  was  operated  in  a  vertical  orientation  spraying  downward  onto  the  heater  substrate. 

Figure  4  shows  a  schematic  diagram  of  the  fluid  loops.  Temperatures  throughout  the  fluid  loops  were 
monitored  using  Type  E  thermocouples  which  were  0.16  cm  (0.0625  in)  in  diameter.  Key  thermocouple  locations 
include  the  entrances  and  exits  of  the  test  chamber,  reservoirs,  and  heat  exchanger  assemblies.  A  positive 
displacement  scavenging  pump  was  used  to  pump  the  FC-72  working  fluid  from  the  sumps  through  a  series  of  three 
air/liquid  heat  exchangers  to  a  reservoir.  Pressure  transducers  were  placed  on  either  side  of  the  pump,  and  the  fluid 
passed  through  a  40  pm  filter  before  reaching  the  reservoir.  The  reservoir  served  as  a  volume  for  the  scavenging 
pump  to  dump  subcooled  FC-72  working  fluid  and  a  pickup  reservoir  for 
the  nozzle  pumps.  Following  the  reservoir,  the  fluid  branched  off  into  two 
identical  flow  paths,  one  for  each  spray  nozzle.  Each  branch  consisted  of  a 
positive  displacement  pump  followed  by  a  15  pm  fdter,  a  flow  meter,  a 
pressure  transducer,  and  a  pressure  switch.  Electrically-actuated  ball  valves 
were  placed  at  the  entrance  and  exit  of  each  branch,  to  permit  isolation  of 
the  test  chamber.  All  testing  was  done  with  subcooled  fluid  at  ambient 
temperature. 

The  experiment  was  operated  and  monitored  with  a  control  panel  and 
data  acquisition  system.  The  data  acquisition  system  consisted  of  a  laptop 
computer  mounted  to  the  top  of  the  experimental  package,  coupled  with  a 
data  acquisition/switch  unit.  Temperature  measurements  were  obtained 


|m  Pressure  transducer 
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Figure  4. 
package. 


Flow  loop  schematic  for  the  experimental  flight 
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Figure  5.  Thermocouple  placement 
within  the  glass  heater  post. 
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throughout  the  chamber,  heater  post,  and  within  the  sump 
annulus  using  sheathed  thermocouples.  Thermocouples  (Type 
E,  0.16  cm  diameter,  time  constant  <  4.0  s)  were  placed  in  the 
chamber  to  measure  chamber  wall  and  vapor  temperatures.  Two 
thermocouples  were  placed  at  the  edge  of  the  heater,  to  measure 
the  liquid  temperature  exiting  off  of  the  heater  surface,  and  in 
the  sump  annulus,  approximately  1.27  cm  below  the  heater 
surface.  Smaller  thermocouples  (Type  E,  0.0254  cm  diameter, 
time  constant  <  0.5  s),  were  mounted  within  the  glass  heater  post 
assembly.  Figure  5  shows  the  geometrical  placement  of  the 
thermocouples  TCI,  TC2,  and  TC3  within  the  glass  heater  post. 
The  dimensionless  radial  and  axial  locations  of  the 
thermocouples  are  given  in  Table  2. 

For  additional  stability  and  accuracy  in  temperature 
measurement,  the  thermocouples  within  the  chamber  and  glass 
heater  post  assembly  were  referenced  to  a  zero  point  reference 
junction  with  a  stability  of  ±0.005°C.  Thermocouples  mounted 
within  the  chamber  and  glass  heater  posts  were  calibrated  using 
an  RTD  temperature  probe  and  a  constant  temperature  bath  in 
the  temperature  range  from  35  to  120°C.  Resulting 
thermocouple  calibrations  yielded  an  uncertainty  of  ±0.02°C. 
Fleater  power  was  monitored  using  voltage  measurements  with 
an  uncertainty  of  ±0.05%  across  the  heater  and  a  precision 
resistor  with  a  resistance  of  0.10  Q  ±  0.02%,  resulting  in  a  heater 
power  measurement  error  of  ±0.07%.  The  flow  rate  of  FC-72 
through  the  nozzle  was  set  at  a  fixed  nominal  flow  rate  and 
monitored  throughout  the  experiment.  The  flow  meters  were 
calibrated,  using  FC-72,  to  an  average  uncertainty  of  ±3.0%  with 
a  maximum  uncertainty  of  ±5. 8%. 


III.  Mathematical  Formulation 


Table  2.  Dimensionless  thermocouple 
locations  within  the  glass  heater  post. 

Thermocouple 
Number 


R 


^post 


TCI 

0.5875 

0.0625 

TC2 

0.5875 

0.1875 

TC3 

0.875 

0.1875 

5 
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Figure  6.  Glass  heater  post  geometry  and 
coordinate  system  used  for  the;  (a)  numerical 
model  and  ( b )  analytical  model  (Not  to  scale). 


With  this  particular  experimental  approach,  the  convective 
heat  transfer  coefficients  around  the  heater  post  assembly  and 
conduction  losses  needed  to  be  determined.  In  order  to  quantify 
the  top  and  wall  convective  heat  transfer  coefficients  and  heat 
lost  due  to  conduction,  numerical  and  analytical  approaches  to 
analyzing  the  heat  transfer  in  the  heater  post  assembly  were 
employed.  Figure  6  shows  the  solution  domain  and  coordinate 
systems  for  the  numerical  and  analytical  analyses.  For  the  most 
general  case,  the  heat  transfer  coefficient  and  heat  flux  on  the  top 
surface  of  the  ITO  heater  may  vary  along  the  radial  and 
azimuthal  directions.  This  is  due  to  the  spray  and  heater  surface  interaction  and  vapor  generation  as  well  as  the  non¬ 
uniformity  of  the  internal  heat  generation  within  the  ITO  heater.  Also,  the  side  wall  heat  transfer  coefficient  could 
vary  along  the  vertical  and  azimuthal  directions.  Both  top  and  side  wall  heat  transfer  coefficients  will  vary  with 
liquid  flow  rate. 

For  the  present  analysis,  the  following  assumptions  were  made: 

1.  The  system  is  under  steady  state  conditions  and  is  axi-symmetric. 

2.  The  thermal  conductivity  of  the  fused  silica  glass  layer,  the  borosilicate  glass  layers,  and  the  ITO  heater 
were  assumed  to  be  identical,  since  the  thermal  conductivities  of  both  glass  types  are  approximately  equal 
at  100°C  ( /ciltl  ~  /Cp()st  ki  ~  A  ~  /c,). 

3.  The  heat  transfer  coefficient  on  the  top  surface  of  the  ITO  heater  is  constant  across  the  radius  where  hMp  (r, 
z  =  IIp,  t)  =  constant. 

4.  The  heat  transfer  coefficient  along  the  outer  radius  of  the  glass  heater  post  assembly  is  constant  where  hm\\ 
(r  =  b,  z,  t)  =  constant. 
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5. 


The  contact  conductances  at  the  interfaces  of  the  glass  layers  are  infinite,  which  means  that  both  the 
temperatures  and  heat  fluxes  at  these  interfaces  are  continuous. 

6.  The  temperatures  used  to  define  the  heat  transfer  coefficients  and  evaluate  fluid  properties  are  obtained 
from  the  experimental  data.  The  freestream  temperature  to  evaluate  the  spray  heat  transfer  coefficient  and 
fluid  properties  is  defined  as  the  average  of  the  fluid  temperature  entering  the  nozzle  and  the  upper  film 
temperature  where  Ta  lop  =  (rspray  +Tuf  )/2  . 

7.  The  freestream  temperature  used  to  evaluate  the  heat  transfer  coefficient  on  the  side  of  the  heater  post 
assembly  is  the  average  of  the  upper  and  lower  film  temperatures  where  Tx  waU  =  (Tuf  +  Tlf  )/2  . 

A.  Numerical  Formulation 

A  steady  state,  axi-symmetric  numerical  model  of  the  cylindrical  glass 
heater  post  assembly  was  analyzed  using  the  finite-element  package 
ANSYS.  The  model  consisted  of  several  layers,  reflecting  the  structure  of 
the  heater,  substrate,  and  post,  as  well  as  thermocouple  locations  within  the 
assembly,  as  shown  in  Fig.  6(a).  Geometric  data  for  the  post  parameters 
shown  in  Fig.  6  are  given  in  Table  3.  The  thermal  conductivity  for  the  entire 
model  was  assumed  to  be  k  =  1.04  W/m-K.  To  reflect  the  symmetry  of  the 
model,  a  zero  heat  flux  boundary  condition  was  imposed  at  the  centerline  of 
the  structure.  The  spray  convection  was  modeled  as  a  convective  heat 
transfer  coefficient,  htop,  on  the  top  surface,  and  the  convection  from  flow 
down  the  side  walls  of  the  post  assembly  was  represented  by  a  separate 
convective  heat  transfer  coefficient  along  the  vertical  surface,  /iwaii.  The  bottom  of  the  post  was  assumed  to  be 
adiabatic.  The  heater  was  modeled  with  a  volumetric  heat  generation  applied  to  the  top-most  layer  (ITO)  of  the 
model. 

For  the  ITO  heater  layer,  assuming  that  the  thermal  properties  are  constant,  the  conservation  of  energy  equation 
is 

^^  +  -8^+8J^  +  -Z-  =  0,0<r<b,L<z<H  (la) 

dr2  r  dr  dz 2  kha  p 

where  g  is  the  internal  heat  generation  per  unit  volume  due  to  ohmic  heating.  The  associated  boundary  conditions 
are 


Table  3.  Geometric  parameters 
of  the  post. _ 


Parameter 

Value 

6 

1.0  pm 

h2 

0.05  ±  0.002  cm 

H 3 

0.1  ±  0.008  cm 

h4 

0.1  ±  0.008  cm 

h5 

9.8  ±  0.05  cm 

L 

10.05  ±0.07  cm 

b 

0.79  ±  0.03  cm 

^^htr  ,  ^wall  rp  ^wall  rp  „  u 

p.  +  ,  1  htr  ,  1  oo, wall  >  r 

or  Htr  Htr 

(lb) 

d^htr  ^top  T  _  htop  _ 

~  '  j  1  htr  ,  1  oo, top  ’  Z  “p’ 

oz  «htr  khtr 

(1c) 

dThtr  ,dT2 

Htr  ~  ~  K2  ~  j  1  htr  ~  1  2  ’  Z  “  Ll- 

oz  oz 

(Id) 

For  the  first  glass  layer,  the  conservation  of  energy  becomes 

d2T?  1  dT ;  d2T? 

+  +  —  0, 
dr2  r  dr  dz2 

(2a) 

with  boundary  conditions 

^  A  h  wal]  T  h  VVa]  |  T  j 

r,  +  ,  2  =  ,  Ao.wall  ,r  =  b, 

or  K.-2  a  2 

(2b) 

u  5rhtr  u  ST 2  rp  rp  _  , 

^htr  n  «  9  1  htr  1 2  ’  z 

OZ  OZ 

(2c) 

dTo  dT \ 

k2^r  =  k3^,T2  =T2,z  =  H2  +ha  +h5. 
oz  oz 

(2d) 
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For  the  second  glass  layer,  the  conservation  of  energy  becomes 


d% 
dr 2 

1  ST ,  d2Tx 

-  +  3  +  3  =0,0<r<b,H4+H5  <z<H3+H4+H5, 

r  dr  dz 1 

(3a) 

with  boundary  conditions 

^^3  ,  ^wall  rp  _  ^wall  7-  „  _  u 

~  '  ,  ^3  ~  j  *  00, wall  ?  '  ~  O’ 

dr  k  ^  k  3 

(3b) 

k2^L  =  h^L,T2  =  T3,z  =  H3+H4+H5, 
dz  dz 

(3c) 

d  T  d  T 

k3-^  =  k4-^,T3=T4,z  =  H4+H5. 

OZ  oz 

(3d) 

For  the  third  glass  layer,  the  conservation  of  energy  becomes 

d  T4  +  ldT4  +d  T4  =0  0<r<b  H  <z<H  +H 
dr2  r  dr  dz2 

(4a) 

with  boundary  conditions 

dT 4  f  ^Wall  rp  _  ^Wall  rp  „  _  u 

~  '  j  *  4  —  ,  *  00, wall  5  '  ~  O’ 

or  k  4  k  4 

(4b) 

dT3  dT4 

*3  ~  -  *4  ~  -  T4,  z  -  H  4  +  H5 , 

dz  dz 

(4c) 

dr4  ST 5 

k4  ,  ~  lc5  ,  14  ~  15’  z  ~  n  5- 

dz  dz 

(4d) 

For  the  fourth  glass  layer,  or 

glass  base,  the  conservation  of  energy  becomes 

^5  +  iSrJ+^J  = 

dr2  r  dr  dz2 

(5a) 

with  boundary  conditions 

^^5  .  ^Wall  rp  ^Wall  rp  7 

..  +  .  T5  ~  j  '00,  wall  ,r-b, 

or  ft  5  « 5 

(5b) 

8T4  dT5 

^4  ~  ~  ,^4  -T5,Z  -  H5, 

dz  dz 

(5c) 

^  =  0,z  =  0. 
dz 

(5d) 

A  grid  independence  study  was  performed  to  ensure  adequate  meshing  of  the  model.  The  coarse  mesh  was 
determined  by  the  fixed  thermocouple  locations.  For  the  coarse  mesh,  the  maximum  element  width  in  the  /'-direction 
was  A r  =  0.00235  m.  In  the  z-direction,  the  heater,  substrate,  and  the  two  thin  layers  were  each  one  element  thick, 
and  the  remaining  portion  of  the  base  was  divided  into  10  elements,  with  a  spacing  ratio  of  30,  so  that  the  element 
size  increased  downward.  This  produced  56  total  elements  in  the  coarse  model.  The  mesh  was  refined  by  doubling 
the  number  of  elements  in  both  directions,  with  the  exception  of  the  heater,  where  its  thickness  remained  at  one 
element.  The  result  was  a  fine  mesh  consisting  of  216  elements.  For  the  grid  independence  study,  typical  convective 
heat  transfer  coefficient  values  of  htop  =  10,000  W/(m-K)  and  /iwan  =  400  W/(m-K)  were  chosen,  along  with  a  heater 
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power  of  Q  =  20  W.  Temperatures  were  recorded  at  all  nodes  and  compared.  The  maximum  percent  difference 
between  the  coarse  and  fine  meshes  was  0.058%. 

B.  Analytical  Formulation 

The  analytical  formulation  was  divided  into  two  regions  as  shown  in  Fig.  6(b).  Region  1  is  defined  as  the  ITO 
heater  with  region  2  being  the  glass  post.  Region  1  was  formulated  as  a  steady-state  one-dimensional  heat 
conduction  problem  with  internal  heat  generation  per  unit  volume  due  to  resistance  heating.  The  top  surface  of  the 
heater  was  treated  as  a  convective  boundary  condition  of  the  third  kind  while  the  lower  surface  was  treated  as  a  heat 
flux  boundary  condition  of  the  second  kind.  Since  the  lower  surface  is  also  the  interface  between  region  1  (ITO 
heater)  and  region  2  (glass  post)  the  interface  condition  must  also  be  satisfied  for  temperature  and  heat  flux.  The 
dimensional  form  of  the  conservation  of  energy  for  region  1  is 


S2Thlr  |  g 

dzhtr  khtr 


=  0,  0  <  z 


htr 


<8, 


(6a) 


with  boundary  conditions 


| 

fehtr 

d^htr 

&htr 


h 

~k 


h 

top  „  _  "top  _ 

J  htr  —  ,  '  co, top 

htr  "htr 

-  ^  t  T  —  T  ~ 
,  ’  1  htr  1  post  ’  z 

"htr 


’  zhtr  “  8 
htr  =  O' 


Defining  non-dimensional  parameters  as 


0  =  - 


T  -  T 


oo,  wall 


1  oo, wall 


the  non-dimensional  form  of  Eqns.  (6a),  (6b),  and  (6c)  become 


— — G  =  0,  0  <  Zht  <  A, 

8ZL 


(6b) 

(6c) 


(7a) 


with  boundary  conditions 


86 

~8Z 


htr 


+  Nu  K0htI  =  Nu  K0x  to p ,  Zhtr  =  A 


htr 

80, 


htr 


8Z 


0htr  ’  ^  htr  ^post  >  -^htr 


htr 


The  non-dimensional  solution  of  Eqns.  (7a),  (7b),  and  (7c)  is  given  by 


#htr  (z)~  ~  —  GZ~  +  tjhV[Z  +  XT  ^tr  +  —  (G8  -  2 tjhtr  )+0x,AoP 


Nu  K 


(7b) 

(7c) 


(8) 


Region  2  was  formulated  as  a  steady  state  two-dimensional  conduction  problem  using  both  a  semi-infinite  rod 
and  a  finite  rod.  The  top  surface  is  the  interface  between  region  1  (ITO  heater)  and  region  2  (the  glass  post) 
satisfying  the  interface  conditions  of  temperature  and  heat  flux.  For  the  case  of  the  finite  rod,  the  lower  surface  was 
treated  as  an  insulated  boundary.  The  side  wall  of  the  heater  post  was  treated  as  a  convective  boundary  condition  of 
the  third  kind. 
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1.  Semi-Infinite  Rod 

For  the  case  when  region  2  is  modeled  as  a  semi-infinite  rod,  the  dimensional  form  of  the  conservation  of 
energy  is 


d  -^post  1  STvost  8  T’post 


dr 2 


-  + — 


r  dr  dz 


=  0,0<r<b,0<  zpost  <  co 


post 


with  boundary  conditions 


drnnct  h 


post  ^  wall  rp  _  ^wall  j, 

'  nnct  rr 


~  j  post  7  J  oo, wall  ’  ^ 

dr  kpost  k„ 


vpost 


dTt 


post  #  \ 


dzr 


’  ^post  ^htr  5  ^post  b. 


(9a) 

(9b) 

(9c) 


"  post  v  post 

Formulating  Eqns.  (9a),  (9b),  and  (9c)  in  the  non-dimensional  form  of  the  conservation  of  energy  results  in 


d2e , 


post  +I^E2!L  +  E$l!L  =  o,  0<  R  <  1, 0  <  Zpost  <  oo, 


8R2  R  8R  8Z 


post 


with  boundary  conditions 


80, 


post 


8R 


+  B impost  -0,  R  -1, 


80, 


post 


8Z 


post 


^post’^post  ^htr’^post 


where  R  =  r/b.  The  non-dimensional  solution  for  the  semi-infinite  rod  becomes 


*postM)=2Bw7postX 


-YvZ 


Jo(rnR) 


i^o(T„)(Bi2+rn2K’ 


(10a) 

(10b) 

(10c) 

(11) 


where  the  yn’s  are  positive  roots  of 


B^o(rn)  =  ra^i(rn)- 


2.  Finite  Rod 

When  region  2  is  modeled  as  a  finite  rod  of  length  L,  the  dimensional  form  of  the  conservation  of  energy  is 


with  boundary  conditions 


d  ^ost 

dr2 


|  1  drpost 

r  dr 


S^post  h 

dr  k 


dT, 


post 


dz 


post 


-  +  ^PL  =  0,0<r<b,0<zpost<L, 

^post 

(12a) 

h 

Wall  rji  _  ' 1  Wall  rji  _ 

post  —  J  co, wall  j  r  ~ 

(12b) 

post  ^post 

Q 1  rji  rji  ___  _  /"v 

7  ’  ^  post  —  htr  ’  ^post  — 

(12c) 

post 
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(12d) 


8T, 


post 


dz 


—  zpost  ~~  L. 


post 


The  non-dimensional  form  of  Eqns.  (12a),  (12b),  and  (12c)  becomes 


d~9  t  1  59 t  8  9 t 

p  -  + - E2!L  + - AA  =  0,0</?<1,0<  ZD0St  <  £ 


57? / 


7?  SR  8Z , 


post 


with  boundary  conditions 


89, 


post 


57? 


+  Bi$Post  —  0,  7?  —  1, 


89, 


post 


az 


post 


^7post  ’  ^post  ^htr  9  -^post 


89, 


post 


az 


=  0,  Zpost  = 


post 


The  solution  to  Eqns.  (13a),  (13b),  and  (13c)  for  region  2  as  a  finite  rod  becomes 


(13a) 


(13b) 

(13c) 

(13d) 


#post(^Z)  =  2Bi,7post  X 

n=l 


where  the  yn’s  are  the  positive  roots  of 


Bi'/o(T„)=5n-/l(Tn) 


(14) 


Analytical  solution  errors  can  result  from  interface  matching  boundary  conditions  (i.e.,  temperature  matching), 
the  number  of  eigenvalues,  and  whether  the  semi-infinite  or  finite  formulation  is  used.  For  this  problem,  the  region  1 
(ITO  heater)  and  region  2  (glass  post)  were  coupled  through  the  interface  conditions  of  heat  flux  and  temperature. 
The  temperature  at  the  top  surface  of  region  2  was  averaged  over  100  evenly  spaced  locations  and  compared  to  the 
temperature  at  the  base  of  region  1.  Temperature  matching  was  performed  by  varying  the  fraction  of  heat  flux  lost 
through  the  base  of  the  heater,  qh  The  criterion  for  matching  the  interface  temperature  was  set  to  less  than  0.0 1°C. 
Solutions  were  calculated  using  30  eigenvalues.  These  solutions  were  verified  to  be  independent  of  the  number  of 
eigenvalues  by  increasing  the  number  of  eigenvalues  to  50.  This  resulted  in  a  maximum  change  in  temperatures  of 
less  than  0.01°C.  In  addition,  solutions  using  the  semi-infinite  and  finite  rod  formulations,  Eqns.  (11)  and  (14),  were 
also  compared  and  found  to  have  differences  of  less  than  0.0 1°C. 


IV.  Results  and  Discussion 


A.  Baseline  Tests 

Terrestrial  (1-g)  baseline  tests  were  run  to  evaluate  and  quantify  the  conduction  losses  through  the  heater  post 
assembly9' 10.  The  bottom  nozzle  with  the  upward  facing  heater  was  operated  using  air-saturated  FC-72  (20-30%  by 
volume)  as  the  working  fluid.  Knowing  the  conduction  losses,  the  heat  extracted  by  the  spray  could  then  be 
calculated  using  a  simple  energy  balance.  Both  the  numerical  and  analytical  formulations  were  used  to  calculate  the 
temperature  distribution  within  the  heater  post  at  three  locations  coinciding  with  the  thermocouple  locations  shown 
in  Fig.  5  and  Table  2.  Numerical  and  analytical  solutions  for  varying  Nusselt  and  Biot  numbers  were  calculated  for 
a  nominal  fixed  flow  rate  (  V  =  5.26><10"6,  8.41><10"6,  and  1 .05x1  O'5  m3/s  (5.0,  8.0,  and  10.0  gph))  and  heat  input  ( Q 
=  20,  30,  and  50  W)  and  then  compared  with  experimental  results. 

Figures  7  and  8  show  examples  of  the  numerical  and  analytical  solutions  compared  to  the  non-dimensional 
experimental  temperatures  91,  01,  and  93,  corresponding  to  TCI,  TC2,  and  TC3,  respectively.  As  would  be 
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Figure  7.  Numerical  results  for  Q  =  50.5  W  at  V  =  Figure  8.  Analytical  results  for  Q  =  50.5  W  at  F  = 
8.41  x  10‘6  m3/s  (8.0  gph)  for  temperature  locations  8.41  x  10‘6  m3/s  (8.0  gph),  for  temperature  locations 
corresponding  to;  (a)  TCI,  ( b )  TC2,  and  (c)  TC3.  corresponding  to;  (a)  TCI,  ( b )  TC2,  and  (c)  TC3. 

expected,  the  temperature  is  least  affected  by  a  variation  of  Biot  number  at  a  location  toward  the  center  of  the  post 
and  closest  to  the  heater,  as  shown  in  Figs.  7(a)  and  8(a),  which  corresponds  to  TCI.  At  locations  farther  from  the 
heater,  i.e.,  TC2  (Figs.  1(h)  and  8(b)),  the  Biot  number  begins  to  influence  the  temperature.  The  greatest  variation  in 
temperature  is  at  the  location  corresponding  to  TC3  (Figs.  7(c)  and  8(c)).  The  temperatures  predicted  by  the 
analytical  solutions  were  consistently  lower  than  those  of  the  numerical  solutions.  This  is  most  evident  by 
comparing  Figs.  7(c)  and  8(c),  where  this  difference  is  the  greatest  with  increasing  Biot  number.  This  is  most 
probably  due  to  how  the  interface  condition  was  treated  in  the  analytical  formulation. 
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was  used  to  determine  the  heat  extracted  by  the  spray.  The  estimated  Nusselt  number  was  determined  from  Fig.  7  at 


the  intersection  of  the  experimental  temperatures,  6\,  02, 
Biot  numbers.  The  local  freestream  temperature,  9X, top, 
was  the  experimentally  measured  film  temperature 
across  the  top  of  the  heater.  The  surface  temperature, 

9  Jr)  was  determined  from  the  numerical  solutions  and 
was  either  treated  as  a  single  point  temperature 
corresponding  to  a  location  directly  above  TCI  or  as  an 
average  across  the  heater  surface.  The  fraction  of  heat 
flux  lost  was  determined  as  the  ratio  of  heat  flux 
extracted  by  the  spray  to  the  input  heat  flux  of  the  ITO 
heater. 

In  a  similar  manner,  the  analytical  solutions  were 
used  to  determine  the  heat  lost  through  the  heater  post 
assembly.  The  estimated  Nusselt  number  was 
determined  from  Fig.  8  at  the  intersection  of  the 
experimental  temperatures,  61,  92,  and  63,  and  the 
analytical  results  for  varying  Nusselt  and  Biot  numbers. 
Since  the  fraction  of  heat  flux  lost,  f  was  an  input 
variable  in  the  analytical  formulation,  the  conduction 
losses  could  be  determined.  Once  a  range  of  Nusselt 
numbers  was  determined,  the  fraction  of  heat  flux  lost 
could  be  estimated  as  shown  in  Fig.  9. 


Figure  9.  Graphical  determination  of  conduction 
losses  for  a  range  of  Nusselt  numbers  for  Q  =  50.5  W 


Figure  10  shows  the  conduction  losses  for  the  three 
approaches.  All  three  approaches  showed  good 
agreement  ranging  from /=  1.0  to  2.5%.  This  result  is 
consistent  with  estimates  of  1.0%  by  Baysinger  et  al.s 
Based  on  these  results,  the  subsequent  analyses  used  a 
fraction  of  heat  loss  due  to  conduction  to  be/=  1.5%. 

B.  Flight  Tests 

Flight  tests  were  conducted  aboard  NASA’s  KC- 
135  Reduced-Gravity  Research  Aircraft  using  air- 
saturated  FC-72  (20-30%  by  volume)  as  the  working 


and  63,  and  the  numerical  results  for  varying  Nusselt  and 
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Figure  10.  Heat  lost  due  to  conduction  through  the  glass 
heater  post  assembly  based  on  the;  (a)  numerical 
solution  based  on  average  surface  temperature,  ( b ) 
numerical  solution  based  on  surface  temperature 
directly  above  TCI,  and  (c)  analytical  solution  (direct 


fluid.  The  bottom  nozzle  with  the  upward  facing  heater  determination  of  f). 
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was  operated  at  nominally  fixed  Weber  numbers  of  We  =  771  ±  19  and  757  ±  15  with  an  inlet  temperature  of  Tspm= 
23.4  ±  0.4°C.  Throughout  the  flight  profile,  the  heater  power  was  varied  while  maintaining  the  test  chamber  at  a 
saturation  temperature  of  Tsi[  =  56  ±  2°C.  Accelerations  generated  by  a  parabolic  aircraft  flight  profile  resulted  in 
the  variation  of  the  non-dimensional  grouping  (Fr1  2Ga)12  =  20  to  66. 

Experimental  data  were  used  to  calculate  the  surface  heat  flux  and  Nusselt  number  using  the  conduction  loss 
through  the  heater  post  of /=  1.5%  as  determined  from  the  baseline  tests.  The  Nusselt  number  was  calculated  using 
the  one-dimensional  energy  balance 
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impost A  2 
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h  (#int  -  # 
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oo, top  ) 


120 


that  accounted  for  the  heat  loss  due  to  conduction  within  the  post. 

Figure  11  shows  the  non-dimensional  heat  input  from  the  heater  surface  versus  the  non-dimensional 
temperature  drop  from  the  heater  surface  to  the  freestream  film.  The  non-dimensional  temperature  drop  increases 
linearly  with  the  non-dimensional  heat 
input  through  ( 1  -f)GA  =  70.  In  this 
portion,  there  was  little  effect  due  to 
varying  (Fr1/2Ga)1/2.  Increasing  the 
heater  input  from  the  heater  surface  to 
( l-f)GA  —110  resulted  in  a  significant 
increase  in  the  non-dimensional 
temperature  difference  between  the 
heater  surface  temperature  and 
freestream  film  temperature  over  and 
above  that  predicted  by  a  linear 
extrapolation.  This  is  most  probably 
due  to  the  initiation  of  vapor  generation 
at  the  heater  surface  and  the  fact  that  the 
freestream  film  temperature  is  used  as  a 
reference  instead  of  the  saturation 
temperature.  The  interface  temperature, 
at  TCI,  for  this  case  was  close  to  the 
saturation  temperature  of  7%  =  56  ± 

2°C.  The  effect  of  (Fr1 2Ga)1/2  is  shown 
by  an  increase  in  the  non-dimensional 
temperature  drop  between  the  heater 
surface  temperature  and  freestream  film 
temperature  with  an  increase  in 


Figure  11.  versus  non-dimensional  temperature  drop  from 

the  heater  surface  to  the  freestream  fluid  for  varying  (Fr1/2Ga)12. 
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19.54,  44.33,  and  65.97. 

The  effect  of  (Fr1/2Ga)1/2  on  the 
Nusselt  number  is  most  noticeable  as 
shown  in  Fig.  12.  The  Nusselt  numbers 
corresponding  to  (Fr1,2Ga)1/2  =  19.54 
and  44.33  were  consistently  higher  than 
the  Nusselt  numbers  corresponding  to 
(Fr1/2Ga)1/2  =  65.97  for  ( \-f)GA  =  30  to 
110.  For  the  non-dimensional  heat 
input  from  the  heater  surface  of  (\-f)GA 
=  110,  the  Nusselt  number,  for  all  of  the 
values  for  (Fr1 2Ga)12,  decreased  most 
probably  due  to  the  initiation  of  vapor  generation, 
continued  to  remain  higher  than  the  Nusselt  number  for  the  (Fr1 2Ga)1/2  =  65.97  case 
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Figure  12.  Nusselt  number  versus  {l-f)GA  for  varying  (Fr1/2Ga)12. 


Flowever,  the  Nusselt  numbers  for  (Fr1/2Ga)1/2  =  19.54  and  44.33 
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Figure  13  shows  Nusselt  number  versus  the  non-dimensional  quantity  (Fr1/2Ga)1/2  for  Weber  numbers  We  =  771 
±  19,  for  (1:/)GA  <  80,  and  757  ±  15,  for  ( l-f)GA  =  107  ±  2.  For  the  experimental  data  shown  in  Fig.  13,  the  Nusselt 
number  tended  to  decrease  with  increasing  (Fr1 2Ga)12. 


V.  Conclusions 

The  effects  of  variable-gravity  on 
single-phase  spray  cooling  with  a  partially 
confined  geometry  have  been 
investigated.  In  general,  the  Nusselt 
number  tended  to  decrease  with  increasing 
(Fr1 2Ga)1/2  for  a  nominally  fixed  Weber 
number.  At  this  time,  the  physics  driving 
the  observed  phenomena  are  not  clearly 
understood.  Some  possible  reasons  may 
be  variation  in  flow  rate,  recirculation  of 
the  spray  film  back  onto  the  heater 
surface,  influence  of  the  confined 
geometry,  effects  of  absorbed  air  in  the 
FC-72  working  fluid,  or  some  unidentified 
thermophysical  phenomena  occurring  at 
the  heater  surface.  Regardless  of  the 
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mechanism,  the  observed  phenomena  appear  to  be  amplified  at  the  onset  of  vapor  generation. 

Based  on  these  very  preliminary  results,  there  appears  to  be  the  potential  for  performance  improvement  for  the 
application  of  spray  cooling  depending  upon  the  magnitude  of  (Fr1 2Ga)1/2.  Flowever,  this  investigation  only 
addressed  the  performance  of  a  single  nozzle  at  a  relatively  fixed  Weber  number.  It  is  not  known  if  the 
implementation  of  a  spray  nozzle  array  will  provide  a  similar  performance  enhancement.  Investigations  with  single¬ 
phase  spray  cooling  for  a  varying  Weber  number  and  the  effect  of  variable-gravity  on  two-phase  spray  cooling  still 
need  to  be  addressed. 
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